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Abstract 



(N 
(N. 

^^ ■ A realistic model for the free rho meson with coupling to two-pion states is 

0^ 



employed to calculate the rho propagator in a hot and dense hadron gas. The 
medium modifications are based on hadronic rescattering processes: interme- 



Oh! 

1^ ■ diate two-pion states are renormalized through interactions with surrounding 



nucleons and deltas, and rho meson scattering is considered off nucleons, 
deltas, pions and kaons. Constraints from gauge invariance as well as the 
full off-shell dynamics of the interactions are accounted for. Within the vec- 
tor dominance model we apply the resulting in-medium rho spectral function 
to compute e"'"e~ production rates from tt^tt^ annihilation. The calculation 
of corresponding e'^e~ spectra as recently measured in central collisions of 
heavy-ions at CERN/SpS energies gives reasonable agreement with the ex- 
perimental data. 
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I. INTRODUCTION 

The main objective of ultrarelativistic heavy-ion coUisions (URHIC's) is to explore the 
phase diagram of strongly interacting hadronic matter, in particular possible transitions 
associated with deconfinement and/or chiral symmetry restoration. Due to it's negligible 
rescattering within the hadronic medium, electromagnetic radiation (photons and dileptons) 
is believed to be the most suitable probe to study the highest excitation states formed in 
the early stages of central collisions. In this context, recent measurements of dilepton in- 
variant mass spectra [|T]-|^ in high-energy reactions of protons and sulfur nuclei with heavy 
target nuclei have raised considerable theoretical discussion. In comparison to the p-A case, 
central S-A collisions show a strong increase in the dilepton yield (after normalizing to the 
total number of charged particles). As several hundreds of pions are produced in the sulfur 
induced reactions, this increase has been attributed to tt+tt^ -^ l^l~ annihilation. However, 
when including this process in numerical transport simulations 0,^] (or hydrodynamical 
approaches |§|,|^) of the collision dynamics, the experimentally observed enhancement in the 
invariant mass range Mi+i- ~ 0.3 — 0.6 GeV/c^ still remains unexplained. 
Several mechanisms that might increase the dilepton yield below the p mass have been pro- 
posed. In particular, the assumption of a a temperature and density dependent dropping 
p mass 1^ as a precursor of the chiral phase transition has been shown to give good agree- 
ment with both the CERES and HELIOS-3 dilepton data @J^. However, before one can 
identify this mechanism as an unambiguous signal of (partial) chiral symmetry restoration, 
the consequences of more 'conventional' scattering processes in the hadronic gas have to be 
well under control. 

In this article we will address this issue, thereby focussing on medium modifications of the 
p meson that are due to phenomenologically rather well established hadronic interactions. 
Such processes typically generate a broadening of the p spectral function and therefore also 
increase the dilepton production at invariant masses below the p mass. 
The coupling of the p meson to intermediate 2-pion states in connection with a modified 



pion propagation in hot/dense matter has been analyzed e.g. in refs. p|-[T3[|. In particular, 
the interaction of the pions with surrounding nucleons and (thermally excited) A's accumu- 
lates substantial strength in the p spectral function for invariant masses below 0.6 GeV/c^. 
Further in-medium scattering processes of the p meson have also been shown to work in the 
same direction, as e.g resonant p-nucleon scattering [|I3| or interactions with pseudoscalar 
mesons 



[LSI . We will here extend on our previous analysis ||T^ by including such processes 
in a consistent off-shell treatment of p propagation in hot hadronic matter. The inclusion of 
the proper off-shell dynamics is of particular importance in the low invariant mass region, 
where the experimentally observed excess of dilepton pairs is most pronounced. 
Our article is organized as follows: in the next sect, we first introduce our model for the 
p meson in free space, which agrees with experimental data on p-wave vrvr scattering and 
the pion electromagnetic form factor in the relevant (timelike) kinematic regime. We then 
briefly review the medium modifications of the p meson generated by a renormalization of 



the pion propagation in a nNA gas as discussed in ref. |TT||T^. In the third sect, we turn to 
the evaluation of in-medium p-baryon scattering. In addition to pN -^ A^(1720), A(1905) 
contributions first discussed by Friman and Pirner |T^, we also account for p-like NN~^, 



AN~^, NA~^, and AA~^ excitations. In sect. 4 we calculate the p selfenergy arising from 
pvr and pK/pK scattering. Within the vector dominance model (VDM) the full in-medium 
p propagator is employed to calculate e"'"e~ production rates in sect. 5. Making use of recent 
transport results to model the temperature/density evolution of central 200 GeV/u S-Au 
and 158 GeV/u Pb-Au collisions, and including experimental acceptance, we finally compare 
our results to the CERES e~^e~ data. 

II. FREE p MESONS AND MEDIUM MODIFICATIONS IN vrvr PROPAGATION 



A. The Model for the p Propagator in the Vacuum 

Our model for p propagator in free space consists of a bare p meson with mass m^"'"'^ 
renormalized by coupling to intermediate 2-pion states. The basic ptttt interaction vertex is 
described by the standard form 

C-p^n = gpnA'^ X d^-K) ■ p^ , (1) 

where tt and p^ denote the isovector pion field and vector-isovector rho meson field, respec- 
tively. After resummation of the tttt bubbles to all orders the scalar part of the p propagator 
can be cast in the form 

Dj(M) = [M^ - {m^-f - T^lAM)]-' , (2) 

which in free space depends on the invariant mass M^ = <lQ~<f' only. The p meson selfenergy 
due to coupling to 2-pion states is calculated as 

S?..(M) = / ^ vp^kf Gl{M, k) , (3) 

with the vacuum 2-pion propagator 

1 1 



'^"^"''^-^^JP^J^WTT,^- "^ = v'"J + *^ W 



and the vertex functions 



v{k) = y o dpniT 2/i; Fp^^{k) . (5) 



The hadronic (dipole) form factor 



, , , 2A2 + m2\2 
F^^(k)= -f—J^) (6) 



accounts for the finite size of the pirn vertex and is normalized to one at the physical reso- 
nance energy mp=0.77 GeV. The subtraction of the p selfenergy at zero energy is necessary 



to ensure the correct normalization ^7^(0)=! for the pion electromagnetic form factor, which 
in the vector dominance model (VDM) is given by 

/',™feare\4 



(M2 - {m^^^^^y - i?eS0(M))2 + (/mS0(M))2 
^ {m'-r \Dl{M)\' . (7) 

Choosing g^^^/4:TT=2.7, Ap=3.1 GeV and m^'^'''^=0.829 GeV yields a satisfactory description 
of the vacuum tttt p-wave scattering phase shifts and the pion electromagnetic form factor 
in the timelike region (cp. fig. 1). 

B. TTTT Propagation in a Hot vrNA Gas 

In general, the p meson propagator in hot/dense matter does not solely depend on 
invariant mass. Since the notion of temperature specifies a certain rest frame of the heat 
bath, Lorentz invariance is broken and the vector meson propagator splits into longitudinal 



and transverse modes (see e.g. ref. Wl 



p Wo, 1)- ^2_ [^harey _ j^L^^^^ q) '^ M^ - {m^-f - Sj(go, q) {wP^'^fM^ ^ ' 
with the standard projection operators 

,u = Oorz/ = 
Pf ={ . . (9) 

5^^ -3^, /x,z/G {1,2,3} 



(the spacelike components of /i and u will be denoted by i and j, respectively). The lon- 
gitudinal and transverse selfenergies are defined by the corresponding decomposition of the 
polarization tensor: 

E^;{qo, q) = S^(go, q) Pi" + Sj(go, q) P^" • (10) 

To calculate medium modifications of the p meson which are generated by dressing the 
intermediate 2-pion states, we adopt the approach of Chanfray and Schuck [jll|, extended to 



finite temperature in ref . |T2| . Restricting ourselves to the case of back-to-back kinematics, 



tlie p meson selfenergy tensor can be written as 

S;i^(go, g = 0) = 5'^ S,..(go, g = 0) . (11) 

Including constraints from gauge invariance the imaginary part of the scalar p selfenergy at 
given temperature T and nucleon-delta density Pna=Pn+Pa takes the form [p]] 

— — ^ Vp^^{kf J [1 + r{ko) + /''(go - ko)] lmD^{qo - /cq, k) 

^ '^ 

* Im{a(go, ko, k)D^{ko, k) + ^^dko, k) + -^T^Tiko, k)}. (12) 

with the longitudinal and transverse spin-isospin response functions 

ni(fco, k) = {kl - {ulf) n°(A;o, k) DL{ko, k) , 

Uriko, k) = {kl - (cu^f) n^iko, k) Driko, k) , (13) 

respectively. The corresponding longitudinal and transverse propagators read 

Ddko, k) = [kl - {tolf - S.(fco, k)]-' , 

DT{ko, k) = [kl - {Lu'^f - C,S.(A;o, k)]-' , (14) 

with Cp=2.2 and 

U''{ko,k) = ^j:^iko,k), (15) 

where J].^{ko, k) denotes the in-medium single-pion selfenergy, i.e. Dl is equal to the single- 
pion propagator D^. The real function a characterizes corrections to the TTTrp- vertex and is 
related to n?j(fco, k) = Ren°(fco, k) by 

a{qo, /Co, k) = l + li^iko, k) + %{qo - ko, k) + ^Ill{ko, k)Ill{qo - ko, k) . (16) 

The bosefactors f^{ko) = [exp(A;o/T) — 1]~^ in eq. (|T2D account for the thermal occupation 
of the pions in the intermediate state. The real part of the p selfenergy is obtained from a 
dispersion integral: 



ReS,..(go) = -Vj -^ , —, . (17) 

^ 

The factor Qq/E'"^ arises from a subtraction at zero energy which ensures gauge invariance 

in the hmit of vanishing three-momentum q. 

The single-pion selfenergy Stt is evaluated within the well-known model of particle-hole 

excitations |18[ extended to finite temperature [|19|. Here the pion interacts with surrounding 

nucleons and thermally excited A's through excitations of the type NN~^, AN~^, NA~'^ 

and AA~^. The corresponding thermal (retarded) Lindhard function in a given excitation 

channel a=ab~^, a,b=N,A, reads 

including the direct (m=l, signe '+') and the exchange (m=2, signe '-') diagram with 

Kk^i^) = i^%A + p'' + k' + 2pkxf/^ . (19) 



For the A width Fa we choose a relativistic parametrization given in Ref. [18|, but neglect 



any in-medium width for nucleons, i.e. Tn = 0. The thermal Fermi distributions, 

r{E;) = (1 + exp[(E; - iia)/T]r' , (20) 

determine the nucleon and A densities at given temperature T and chemical potential /iat, 
yUA to be 

P^(T)=4/^r(Ef;/i^,T) 

PAm = 16/^/^(i?f;/^A,T). (21) 

Throughout this article we will assume chemical equilibrium characterized by a common 

baryon chemical potential /iB=/iAr=/iA and /i7r=0. 

From the Lindhard functions one obtains the so-called pionic susceptibilities Xa^ as p 



xL°)(^,fc) = (l^^IlA^y siia) Uu^,k) , (22) 

with spin-isospin factors SI (a), a standard monopole form factor 

(A^=1.2 GeV) and coupling constants f-„a related via fnNA='^fTrNN (Chew-Low factor pU[ ) 
and f-K/^/\=\f-KNN (constituent quark model estimate 0]), cp. table |. For a realistic 
description of the pion selfenergy the bare susceptibilities have to be corrected for short- 
range correlation effects between particle and hole. These are parametrized in terms of 
Migdal parameters g' api leading to a system of coupled equations 

x. = xi°^ - E xL°) g\, X, . (24) 

/3 



It is solved by an elementary matrix inversion with the final result [|19| 

T.^{uj,k) = -k'' Y.Xa{^,k) . (25) 

a 

The Migdal parameters are fixed to g' ^^^ = 0.8 for af3=aa~^hh~^ and g' ^f^ = 0.5 for all others. 
With the diagonal form eq. ( [TTl) of the p selfenergy tensor the corresponding in-medium 
propagator can be written as 

^;'(go, g = 0) = 5'^ D,{qo, q = 0) , (26) 

where, in analogy to the free case eq. @, the scalar part is simply 

D,{qo, q = 0) = W - {m';n' - ^p^AQo, q = 0)]-' . (27) 

In fig. 2, real and imaginary part of Dp are displayed for different temperatures T at a fixed 
baryon chemical potential of yUAr=/iA=0.39 GeV (this value is chosen with respect to what 
will be used for calculating dilepton spectra in central S-Au collisions in sect. 5). With 
increasing temperature (which at fixed ^n, Pa corresponds to a simultaneous increase in 
Pna=Pn+Pa) we find a strong broadening of the p spectral function. Although the peak 
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is moderately shifted to higher energies (in accordance with the findings of refs. P-p!T[|), 
an appreciable enhancement over the vacuum result is observed for invariant masses below 
M=0.6 GeV. 

In the calculation of e^e~ production rates from in-medium tt+tt" annihilation (sect. 5) 
the p spectral function will be weighted with an additional factor of 1/M^ steming from 
an intermediate photon propagator. It is instructive to see how this combines with our 
results for Im Dp (cp. fig. 3). At highest temperatures the major part of the strength 
resides in the invariant mass range M=0. 1-0.5 GeV, which simply refiects the softening of 
the pion dispersion relation: the low-lying NN"^ and AA^^ modes lead to a population 
of invariant masses around M=0.2 GeV, whereas the AA^^^ excitations, strongly washed 
out by thermal motion and the large A width, show up at somewhat higher M~0.4 GeV. 
In addition, the shift of strength to low invariant masses is accelerated by the appearance 



of the pion bosefactors in eq. (0). Even at lower densities one clearly recognizes various 
combinations of pionic branches in the 1/M^- weighted p spectral function. 
To investigate the role that the pionic collectivity plays for the build-up of the medium 
effects just discussed, we replaced the full propagators Dt,=Dl and Dt, entering eq. (P^, 
by their vacuum forms and recalculated the corresponding p propagator, cp. fig. 4: without 
the pion collectivity the broadening of the p spectral function is strongly reduced, and the 
peak is now slightly shifted downwards (short-dashed line). For small temperatures this 
trend becomes even more pronounced (dotted line). A similar behavior has recently been 



found in a calculation to lowest order in the nuclear density [^ 



III. RHO SCATTERING IN BARYONIC MATTER 

Besides the medium modifications caused by dressing the intermediate 2-pion states, di- 
rect interactions of the p meson with surrounding hadrons in the gas have to be considered. 
In analogy to the in-medium behavior of pions one might expect strong effects from p scat- 
tering off baryons (provided their abundance in the hadronic gas is appreciable as seems the 



case in central URHIC's at CERN energies p3[)- ^^oia meson exchange models such as the 
Bonn potential |^ one certainly knows that e.g. the pNN (or pNA) coupling constant is 
quite sizeable, even though the corresponding s-channel process pN —>■ N —>■ pN is kinemat- 
ically strongly disfavored. Such a kinematic suppression will be much less pronounced with 
increasing energy of the resonance in the intermediate state. Indeed, there are at least two 
well established resonances in the particle data table |2^ which strongly couple to the pN 
decay channel, namely A^(1720) with isospin/spin-parity /(J^)=|(| ) and a branching ra- 
tio rjv(1720)-piv/r*^*(i720) > 0.7, and A(1905) with /(J^) = |(| + ), rA(1905)-piv/r*4i905) > 0.6. 

This led Friman and Pirner to the idea [Q of considering p-like particle-hole excitations 
('rhosobars') in nuclear matter of the type pA^(1720)A^~^ and pA(1905)A^~^. In analogy to 
the well-known 'pisobar' (ttAA^"^), the (non-relativistic) interaction Lagrangian for the two 
'rhosobars' can be written as 

nip 2 

Cp^'N = ^^^^^ ^l SI' q' (T ■ pi) ^M + h.c. (28) 

TTlp 2 

{N* = Ar(1720), A* = A(1905)). Here, Ss {S7) denote the spin ^ ^ f (f) transition 

2 2 

operator (tensor), r, T the isospin | and the isospin | — * | transition operator, respectively, 
and g the 3-momentum of the p meson. The coupling constants fpB*N {B*=N* or A*) are 
adjusted to reproduce the partial decay widths Tb*^pn- To obtain realistic values for them 
it is important to account for the finite width of p meson, since both resonances (especially 
the A^(1720)) lie rather close to the pN threshold (the assumption of a sharp p meson would 
allow for only a small phase space, which in turn would result in an overestimate of the 
coupling constant). Starting from the general expression for the differential two-body decay 



width, given e.g. in ref. [^, we arrive at 



where SI{N* -^ pN)=2, SI{A* -^ pN)=l/3 are initial-state averaged, final-state summed 
spin-isospin factors and FpB*N{.<icm) is a hadronic monopole form factor {KpB*N='^ GeV) 

10 



depending on the p/N decay 3-nionientum in the resonance rest frame, 

Hem — T~ • K'^^) 

The p meson spectral function, 

Ap{M) = -2 lmDp{M) , (31) 

is taken from eq. (0) and the upper integration hmit is Mmax=\'^ — mj\f. With a branch- 
ing ratio of 70% for the A^(1720) decay into pN and a total width of 0.15 GeV at the 
resonance energy m^r* =1.72 GeV, eq. ( pQ]) determines the coupling to be fpj^*j^/4:7!'=6.99. 
Applying the same procedure for the A(1905) with a branching ratio to pN of 60% we find 
/pA*Jv/4vr=10.64. 

Having fixed the coupling constant we can now compute the in-medium selfenergy tensors 
for p-like B*N~^ excitations from the interaction vertices of eq. (^), leading to 

* E (^A^..|(S| X q)^\h^) {\\^\{S\ X gTI^A^v*) (32) 

and a similar expression for the pA*N^^ bubble involving the 5*5 tensor. Both the isospin 
factors 

IF{pN*N-^) = 2 

IF{pA*N~^) = - (33) 

o 

and the Lindhard functions 

(a = A^*A^~^, A*A^~-'^) are equivalent to the pionic case, eq. (pISl). The total widths are each 
taken as the sum of pN channel and vrA^ channel, 

Tfi. = '^B*^pN + Ts'^TrAr , (35) 

11 



where the ttB*N couphng constant is chosen such that the total width matches it's experi- 
mental value at the resonance mass. 

The transversality of the pN*N coupling automatically ensures that the vector current is 
conserved, i.e. 

q, E^y^., = , (36) 

which, when coupling to the photon, is a necessary condition for gauge invariance. Thus, for 
the A(1905) contribution we only take into account the transverse part. The corresponding 
scalar parts of the selfenergy can be extracted by applying the transverse projection operator 
from eq. @: 



J- / J pa J^pa{Q) \ ciT( \ 2 



2 y rap 



2 

SI{pa) q^ (ppaiqo,q) 



9'xS(?o,g), (37) 



where the spin-isospin factors are summarized in table ||. Guided by the experience one has 
from pion nuclear physics, we also account for short-range correlation effects in the p-like 
particle-hole bubbles, parametrized by Migdal parameters g'. The renormalized selfenergy 
in a given excitation channel then becomes 



Sp.(go, q) = -q' . , ;" w'r . • (38) 

For more than one channel one again has to solve a matrix problem as given by eq. 



Xpa = Xfl - E Xfl g'ap Xp0 . (39) 

Both the A^(1720) and A(1905) resonances have sufficiently high masses which allow to 
disentangle their coupling strength to the pN channel through the corresponding decay 
process. However, within our full off-shell treatment of the in-medium p propagator we 
are also able to incorporate lower-lying p-like particle-hole excitations. The most obvious 
candidate is, of course, the A(1232): it's typical excitation energy of about 300 MeV closely 

12 



coincides with the invariant mass region where the experimentally observed excess of dilepton 
radiation in heavy-ion collisions at CERN energies sets in |1|,^ . An appropriate Lagrangian 
is given by 

Cp^N = ^ ^i iS._ X g) ■ (T ■ p) vl/^ + h.c. , (40) 

which automatically includes pNA~^ excitations as well. To complete the picture of low- 
lying p-baryon excitations we furthermore account for NN~'^ and AA^^ bubbles described 
by 

rup 

CpAA = ^^^^i (2S3 X g) ■ (6 • p) M/A , (41) 

nip 2 

where Sa and G denote the spin-| and isospin-| operators, respectively. The coupling con- 
stant fpNN is taken from the Bonn potential, and fpNA='2fpNN, fpAA=lfpNN, cp. table 0. 
Including Migdal parameters as in the pionic case (all additional ones are uniformly set to 
g'^ij=0.5), the total selfenergy for p-baryon interactions in the particle-hole picture becomes 

^pBB~^{qo,q) = -q^ J2Xpa{qo,q) (42) 

a 

where the summation is over a = NN''^, AN''^, NA''^, A A'^, N (1720) N^^^, A{1905)N-^, 
and the full susceptibilities are determined from solving eq. (|39|). 
In fig. 5 we display the transverse part of the p spectral function, 

^"'"' '^' '' ^ '"" ( M^-«^.p-E/(M)-W.te.J • («' 

at fixed baryon chemical potential and with no medium effects included in the 2-pion states. 
With both increasing density (upper panel in fig. 5) and 3-momentum (lower panel in fig. 5) 
the 'rhosobar' branches build up an appreciable enhancement for invariant masses below 
M ~0.6 GeV, together with a depletion of the p peak. However, due to the thermal motion 
of the nucleons and the relatively large widths of the A^(1720), A(1905) and A(1232) res- 
onances, the overlapping branch structures are hardly visible; only at the highest nucleon 
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density considered (piv=0.7po at T=0.17 GeV, dotted line in the upper panel of fig. 5) one 
recognizes a shoulder around M=0.6 GeV corresponding to the piV*iV-channel. 

IV. RHO SCATTERING IN A THERMAL MESON GAS 

In central collisions of heavy nuclei at CERN energies (160-200 GeV/u) large numbers 
of secondary mesons are produced, clearly exceeding the number of primary nucleons. How- 
ever, as one knows from the pion, it's properties are typically much less affected in a hot 
meson gas than in the nuclear medium at comparable densities 1^^. Nevertheless, a reli- 
able description of p meson properties at CERN-SpS energies certainly has to account for 
interactions with surrounding mesons. An analysis of collision rates for (on-shell) p-mesons 
in a tt/K gas has been performed in ref. [jl5[, where a moderate collisional broadening of 



r;;°"(T<0.17 GeV)<60 MeV has been found. 

Our point here is to include pvr and pK scattering as an additional contribution to the (off- 
shell) p selfenergy in hot hadronic matter. In analogy to the previous section we assume the 
dominant contribution to arise from s-channel resonance formation, i.e. pn —>■ ai(1260) and 
pK/pK -^ i^i(1270)/-ft'i(1270). For the corresponding interaction Lagrangian we choose 
the one proposed in ref. [^, which is (isospin structure suppressed) 



-^Trpai = ^Trpai Ct^ id'"' Q ' Ptt ' ?>^) Pu TT , (44) 

being compatible with gauge invariance within the VDM approach, a^, pu and vr denote the 
fields for the ai(1260) meson, rho meson and pion, and q, p-,^ the rho and pion 4-momenta, 
respectively. Similarly, we describe the K pKi{127Qi) interaction by 

CkpK, = GkpK, b^ {g^^ q ■ Pk - q^Px) Pu K (45) 

{K: kaon field, b^: 7^1(1270) field). Along the same lines as in the previous section we fix 
the coupling constants G by adjusting the partial decay widths Fa^^^p, Tki^kp to their 
experimental values. Including the finite width of the p meson one has 
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GL„. ^'TMdM 



^».^-'(^)^S7/^-^<-) 



2m 

* q,^ [^{s - M2 - miy + M\ml + qlJ] F^p^Mcmf (46) 

with a hadronic dipole form factor FT,pa^{qcm) {^■Kpai='^ GeV) and Mmax = y/s — m.„. As- 
suming Ta,-,np{^/^ = ma,) = C/ =0.4 GeV with m„^ = 1.23 GeV yields G^pa, = 18.7 GeV-S 
which is not very different from the value of 14.8 GeV~^ obtained for a zero-width p meson in 



ref. [^. An analogous calculation for the Ki width, taking TKi^Kp{y/s = mKi)=O.OQ GeV 



and mK, = l-27 GeV @, gives Gkpk^=15.3 GeV~\ 



Within the imaginary time (Matsubara) approach the p meson selfenergy tensor arising 
from interactions with surrounding pions at temperature T can be calculated as 

Sp.^(^o, q}=J (03 ^[/^K) - r'K + ^o)] M^;ip^, q) (47) 

with bosefactors f", f^^ and the isospin averaged forward scattering amplitude M^^. As- 
suming the latter to be dominated by ai resonance formation one derives from eq. (^41): 



(2/ai + 1) 2 p ( \2 



M^;,^ip^,q) = )^;'^^^ Gip,^ F^ 

<9''' P.-q- P'.q') Da.Mk) {g'^'q " P. - qV.) , (48) 



where the Oi propagator at 4-momentum k = (p^ + q) is given by 



^'^i.'^a(^) = ^^''\ '^. '^ ^ll . 49 



(s = kxk ). The total ai width r„° (s) is taken from eq. (^61) . Using the projection operators 
defined in eq. @ we can extract the explicit form of the longitudinal and transverse part of 
the selfenergy as 

Sp^.a,(?o,g) = (PL);.. s;^l,(go,g) 



*- ^ '1 "1 



s-m\ +im«,r^f 
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^l^Ml) = Upt),u ^'pLAlo,<l) 



I ^-- / (i^ f^^^<) - ^"^"-' + ^°^] ^--^^^ 



[l(.-M2-m^)2^-p-^M2(l-^)(l-a:2)] 



(50) 



with X = cos 9, = l{p,q). Completely analogous expressions are obtained for Kp scattering 
via ^1(1270) resonance formation (replacing tt ^ K, ai ^ Ki, note that the isospin 
averaging factor in eq. (|48| ) changes from 1 to 2/3 and that the total width entering Dk^ is 
taken to be T^^^_^=2Tki^kp'i for antikaons one has T,pj^f^^=T,pKKi as long as Px=iik)- The 
spin-averaged p propagator in a hot t^-K gas is then given by 



l[-9,u + ^-0^]D^'^{M,q) 



1 



3 VM2 - (m^-)2 - SO,,(M) - S^.„,(M, q) - 2S^^^^(M, g) 

+ M2 _ (m^arey _ ^o^^^M) - S J,„, (M, g) - 2Sj^^^(M,g) j ' ^^^^ 

the imaginary part of which is shown in fig. 6 for /i^r = px = 0. With increasing temperature 
(upper panel) np and Kp/Kp scattering simply lead to a moderate broadening of the spectral 
function, which is similar to the results of ref . [^ . The impact of a finite 3- momentum is 
rather small (lower panel of fig. 6). 

V. e+e SPECTRA IN CENTRAL HEAVY-ION COLLISIONS AT CERN/SPS 

ENERGIES 

The general expression for the dilepton production rate 'R=dNi+i- / d*x per unit of 4- 
momentum d!^q in a hadronic medium of temperature T can be written as 

|^ = L,.fe)i^-(,). (52) 

To lowest order in the electromagnetic coupling a=l/137 the lepton tensor is given by 
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where we neglected the rest mass of the leptons compared to their individual 3-momenta 
|p+|, |p„| (in the following we will focus on the e+e^ case), q = P+ + P- is the total 4- 
momentum of the pair in the heat bath. 

To calculate the dilepton prodcution rate from the tt+tt^ -^ e^e^ annihilation process we 
employ the phenomenologically well established vector dominance model (VDM). It relates 
the hadronic part of the electromagnetic current to the third component of the isovector p 
meson field as 

J' = ^^\^P3 , (54) 

with the universal VDM coupling constant g = Qp-mv As a consequence, the hadronic 
tensor H^^ can be expressed in terms of the imaginary part of the retarded p propagator in 
hot/dense matter: 

m^{q,, q; ps, T) = -/^ (go; T) ^^\^ 11X1^^(^0, g1 /iB, T) (55) 

{f^(qo,T) = (e'^"/"^ — 1)~^)- Oil account of current conservation, qfjliaD^'^ = 0, and using 
the decomposition eq. (§) we arrive at 

^A^^Aq = ^^^ Jp- 3 (I"^^P (^0, q; /XB, T) + 2 ImD^ (go, g; Pb, T)) (56) 

with 

T ni,r. ^^ ImS^'^(go,g;/iB,r) 

ImD (go, g; /x^, T) = — J -jj^ — - . (57) 



|M2-(m^--)2-S^'^(go,g;/iB,T)| 



2 



The full p meson selfenergy is obtained from combining the results of the three preceeding 
sections, eqs. (|l|), (0), (||) and (11): 

Sp = T^p-mr + ^pT^a^ + '^pKKi + '^pKKi 

Sp = Sp^^ + T^pBB-^ + Sp^^^ + Sp^^^ + '^pKKi ■ (58) 
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In fig. 7 we show our final result for the spin averaged p spectral function Im Dp = |(IniD^ + 
2IniZ)J), which does not exhibit any principly new aspects as compared to the analyses in 
the previous sections. 

For calculating e^e~ invariant mass spectra as measured in the CERES experiment the 
differential rate eq. ( |56D has to be integrated over 3-momentum as well as the space-time 
history of a central 200 GeV/u S-Au reaction. Therefore we make use of recent transport 
calculations [^] according to which the corresponding temperature evolution can be well 
represented by an exponential decrease 

T(t) = (r - T°°) e"*/^ + T°° (59) 

with an initial temperature T*=0.170 GeV, r=8 fm/c and T°°=0.110 GeV. The initial 
baryon density is predicted by RQMD simulations to be p^^2.5po [ID- When including 
the 12 lowest lying baryonic resonances this corresponds to a baryon chemical potential of 
/iB=0.39 GeV, which is slightly lower than the value we employed in |]12[. If we furthermore 



assume chemical equilibrium, i.e. fiB=const and fimeson=0, the time evolution of baryon and 
meson densities is completely determined by eq. (pQ]). The such obtained baryon density 
psit), e.g., turns out to be in good agreement with the transport simulations. After a 
variable transformation qodq^ = MdM and with an isotropic density profile of the hadronic 
fireball at each stage of the expansion, eq. ( ^6|) becomes 



* fd^q — nqo;nt)) lmDp{M,q;pB,T{t)) A{M,q) . (60) 

A{M, q) accounts for the experimental acceptance cuts on the dilepton tracks as applied in 
the CERES detector {i.e. pt>0.2 GeV, 2.1<r7<2.65 and ee+e->35 mrad). We here also 
include it's finite mass resolution by folding our theoretically calculated spectrum with a 
gaussian 

CM, m = ^= exp - \ , ^ J^ , (61) 



where the width a{M) is extracted from ref. ^. The dimensionless normahzation constant 
A'"o is fixed using free p decays by matching our corresponding total (M-integrated) e~^e~- 
yield with the transport resuhs of Li, Ko and Brown 0. In eq. ( |60D t/o denotes the freezeout 
time of the hadronic firebaU. For central S-Au reactions we use tfo=10 fm/c (corresponding 
to a baryonic freezeout density pB(tjo)=0.32po), but our results are rather insensitive with 
respect to moderate variations in tfo. The fireball volume V{t) is determined by the baryon 
density as 

where Nb denotes the number of participating baryons (which can be absorbed in the 
normalization constant A^o)- 

In fig. 8 we compare the resulting dielectron spectrum, supplemented with contributions from 
free Dalitz decays {'n'o,V ~^ 7e"'"e~, u -^ Tr^e^e^) and free u -^ e~^e~ decays as extracted 
from ref. |Q , with the CERES data from central S-Au collisions at 200 GeV/u. Apparently, 
the combined medium effects in the p propagation lead to quite reasonable agreement with 
the experimental spectrum. As was identified in the preceeding sections, both the dressing 
of the intermediate 2-pion states as well as particle-hole-type excitations by the p meson 
('rhosobars') constitute the major part of the effect. 

We also performed a similar analysis for central Pb-Au collisions at 158 GeV/u. Here, 
RQMD predicts initial temperatures and baryon densities of T*^0.180 GeV and p^^4po |p3|| , 
which corresponds to a baryon chemical potential of yUB=0.408 GeV. Taking r=10 fm/c and 
T°°=0.105 GeV in eq. (p9D, the time evolution of temperature and density as found in 
transport calculations [|29| can again be well reproduced. The much larger system size as 
compared to the S-Au case increases the fireball lifetime substantially. We account for this 
by choosing a freezeout time of tfo=20 fm/c corresponding to a freezeout temperature and 
density of T'^°=0.115 GeV and pB{tfo)=0.18po, respectively. This is somewhat lower than 
in S-Au collisions, which refiects the fact that the hadronic cooling becomes more efficient 
with increasing system size. The comparison of the resulting dielectron spectrum with 
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preliminary data from the CERES/NA45 collaboration confirms our findings for S-Au case: 
hadronic rescattering in in-medium p propagation seems to resolve the discrepancy between 
the experimental data and theoretical results based on the free tt+tt^ annihilation process. 

VI. SUMMARY AND CONCLUSIONS 

Starting from a realistic model for the p meson in free space we have studied p propagation 
in hot hadronic matter and it's impact on dilepton production in URHlC's at CERN-SpS 
energies. Since the free p strongly couples to vrvr states, an important medium effect arises 
from the modification of the pion propagation: a strong softening of the pion dispersion 
relation, generated by the coupling to AA^^^ as well as low-lying NN^^, AA^^ excitations 
and further weighted by thermal occupation factors in the hot gas, leads to an apprecia- 
ble enhancement of the p spectral function from zero to about 0.6 GeV in invariant mass. 
However, as we have shown earlier [|1^], these modifications do not yet allow for a quanti- 



tative description of the CERES data. We furthermore included direct interactions of the 
p within the hot hadronic enviroment. Corresponding scattering processes off surrounding 
nucleons, deltas, pions and kaons, were assumed to be dominated by s-channel resonance 
graphs, which we derived in accordance with constraints from gauge invariance. Within 
a full off-shell treatment the most notable contributions come from 'rhosobar'-type AA^~^, 
A^(1720)A^~^ and A(1905)A^^^ excitations, which result in a further shift of strength to low 
invariant masses in the p spectral function. ai(1260) and 7^1(1270) resonance formation in 
pvr and pK scattering turned out to be less significant. Within the VDM we applied our 
full p propagator to compute e~^e~ production rates from in-medium vrvr annihilation. These 
were integrated according to a temperature/density evolution of central S-Au (200 GeV/u) 
and Pb-Au (158 GeV/u) collisions as found in recent transport calculations, including exper- 
imental acceptance conditions. It turns out that the observed enhancement in the CERES 
data can be accounted for. 
One can certainly think of other medium effects in p propagation that are not included in 
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our present study. Among those are, e.g., finite temperature modifications of the VDM {i.e. 
corrections to the pj vertex), or a-tadpole graphs [^, which represent the coupling of the 
p meson to a scalar field generated by the surrounding hadrons. Furthermore, the vector 
dominance assumption for the nucleonic sector is presumably less accurate than for the pion, 
which might overestimate the effects of the baryonic medium somewhat (especially for very 
low invariant masses). However, we do not expect these shortcomings to affect our results in 
a major, qualitative way. Thus, it seems justified to conclude that medium modifications of 
the p meson arising from hadronic rescattering can essentially explain the low-mass dilepton 
enhancement as found in recent experiments at the CERN-SpS. After all, the Brown-Rho 
conjecture of a dropping p mass might in fact be a related phenomenon, i.e. in some sense 
an efficient way to parametrize the shift of strength to low invariant masses in the p spectral 
function. However, in our dynamical calculation this shift of strength is due to a strong 
broadening of the resonance, whereas BR-Scaling predicts a downwards shift of the peak 
caused by a strong scalar mean field in the hadronic medium. Some of this discrepancy 
might be resolved when including the a tadpole graph in our framework. 
Further refinements also require to include a finite g-dependence of the 2-pion selfenergy. 
Finally, we need to improve our treatment of the URHIC-dynamics; here, a hydrodynamic 
approach, which accounts for the full off-shell behavior implicit in our dilepton production 
rate, is mandatory. Work in all these directions is in progress. 
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TABLES 
TABLE L Spin-isospin transition factors and coupling constants for pion induced (longitudinal) 

p-wave particle-hole excitations in a hot NA gas. 

na ttNN-^ ttAN-'^ ttNA-^ ttAA''^ 

SI{iTa) 4 16/9 16/9 400 

/2^/47r 0.081 0.324 0.324 0.00324 

TABLE IL Spin-isospin transition factors and coupling constants for p meson induced (trans- 
verse) p-wave particle-hole excitations in a hot NA gas. 

pa pN{1720)N-^ pA{l905)N-^ pAN'^ pNA'^ pNN'^ pAA'^ 

SI{pa) 16/3 8/5 32/9 32/9 8 800 

/2„/47r 6.99 10.64 18.72 18.72 4.68 0.19 
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Figure Captions 



Figure 1: Our fit to the p-wave vrvr scattering phase shifts (upper panel) and the pion electromag- 
netic form factor in free space (lower panel); the squares in the lower panel correspond 
to the Gounaris-Sakurai formula [0, which gives an accurate description of the data. 

Figure 2: Impact of in- medium vrvr states on the p propagator (upper panel: imaginary part, 
lower panel: real part) at zero 3-momentum in a hot irNA with chemical potentials 
fiB=fJ'N=fJ'A=0-'i^ GeV and fin=0- The chosen temperatures of T=0.170 GeV (dotted 
lines), T=0.149 GeV (short-dashed lines) and T=0.127 GeV (long-dashed lines) cor- 
respond to A^A densities of pna/po=^-4:1, 0.63 and 0.32, respectively. The solid lines 
represent the results in free space. 

Figure 3: Imaginary part of the p propagator, weighted with a photon propagator 1/M^, in a 
hot ttNA gas; medium modifications are due to the dressing of the intermediate 2-pion 
states; line identification as in fig. 2. 

Figure 4: Imaginary part of the p propagator with in- medium vrvr states at normal nuclear mat- 
ter density; 

long-dashed line: full result at T=0.150 GeV; short-dashed line: without pionic col- 
lectivity (as described in the text) at T=0.150 GeV; dotted line: without pionic col- 
lectivity at T=0.050 GeV; full line: free space. 

Figure 5: Imaginary part of the in-medium p propagator when accounting for pN and pA scat- 
tering in a hot baryon gas at chemical potential pb=Pn=Pa=0.39 GeV; 
upper panel: for fixed 3-momentum g=0.5 GeV and temperatures T=0.127 GeV (long- 
dashed line), T=0.149 GeV (short-dashed line) and T=0.170 GeV (dotted line); 
lower panel: for fixed temperature T=0.149 GeV and 3-momenta g=0.25 GeV (long- 
dashed line), g=0.5 GeV (short-dashed line) and g=0.75 GeV (dotted line). 



26 



Figure 6: Imaginary part of tlie in-medium p propagator when accounting for pvr and pK/ K 
scattering in a hot meson gas at zero chemical potential; 

upper panel: for fixed 3-momentum g=0.5 GeV and temperatures T=0.127 GeV (long- 
dashed line), T=0.149 GeV (short-dashed line) and T=0.170 GeV (dotted line); 
lower panel: for fixed temperature T=0.149 GeV and 3- momenta g=0 (long-dashed 
line), g=0.5 GeV (short-dashed line) and g=l GeV (dotted line). 

Figure 7: Our final results for the imaginary part of the in-medium p propagator after inclusion 
of in-medium vrvr states as well as p meson scattering off nucleons, delta's, pions and 
(anti-) kaons. The curves correspond to fixed 3-momentum (g=0.5 GeV) and chemical 
potentials (/iB=0.39 GeV, Pmeson=^) at temperatures of T=0.127 GeV (long-dashed), 
T=0.149 GeV (short-dashed) and T=0.170 GeV (dotted). 

Figure 8: Dielectron spectrum from central S-Au (200 GeV/u) collisions; the dots are the 
CERES/NA45 data with statistical errors (vertical bars) and systematic uncertainty 
added independently (resulting in the thick horizontal bars); the dotted curve is the 
sum of free Dalitz (dashed-dotted line) and uo decays (both extracted from recent 
transport results |^) as well as free it~^it~ annihilation; the full curve is obtained when 
evaluating the tt+tt^ contribution with the medium modified p propagator. 

Figure 9: Same as fig. 8, but for central Pb-Au (158 GeV/u) collisions, where the dots are 



preliminary data from the CERES/NA45 collaboration ||3y]; here, only the statistical 
errors are displayed. 
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